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Exit probability in a one-dimensional nonlinear q-voter model. 
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We formulate and investigate the nonlinear g-voter model (which as special cases includes the 
linear voter and the Sznajd model) on a one dimensional lattice. We derive analytical formula for 
the exit probability and show that it agrees perfectly with Monte Carlo simulations. The puzzle, 
that we deal with here, may be contained in a simple question: "Why the mean field approach gives 
the exact formula for the exit probability in the one-dimensional nonlinear g- voter model?". To 
answer this question we test several hypothesis proposed recently for the Sznajd model, including 
the finite size effects, the influence of the range of interactions and the importance of the initial step 
of the evolution. On the one hand, our work is part of a trend of the current debate on the form 
of the exit probability in the one-dimensional Sznajd model but on the other hand, it concerns the 
much broader problem of nonlinear g-voter model. 

PACS numbers: 64.60. Ht, 05.70.Ln, 05.50.q 
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I. INTRODUCTION 

Linear voter model [1| , one of the most recognized in 
a field of non-equilibrium phase transitions, is not only 
a toy model of an Ising spin system but also caricature 
of opinion dynamics. One of the main reasons of its im- 
portance is the fact that the linear voter model (VM) is 
solvable in arbitrary spacial dimension. However, from 
social point of view it is definitely to simplified and there- 
fore several other models of opinion dynamics, based on 
Ising spins, have been introduced (for excellent resent re- 
view see [2]), e.g. Sznajd model Q or majority model 
[J, Q . Of course, it happens that a seemingly different 
models give the same results or even can be formulated 
in such a way that appear to be identical. For example, 
it has been shown that the original one-dimensional Sz- 
najd model [3| can be rewrite as a classical voter model 
[^ . However, the most commonly used version, in which 
only the unanimous pair changes the state of the system, 
differs significantly from VM. Nevertheless, it has been 
suggested that this case could be described by a broader 
class of nonlinear voter models 7, 8]. 

Recently, particularly interesting nonlinear variant of 
the voter model, the g-voter model has been introduced 
[8|. In the proposed model q, randomly picked (with pos- 
sible repetitions), neighbors influence a voter to change 
opinion. If all q neighbors agree, the voter takes their 
opinion; if they do not have a unanimous opinion, still a 
voter can flip with probability e. For q = 2 and e = the 
model is almost identical with Sznajd model on a com- 
plete graph [9| . The only difference is that in the g- voter 
repetitions in choosing neighbors are possible. However, 
for q — 2 and reasonable large lattice size this difference 
is negligible. In this paper we formulate and investigate 
the g- voter model on a one dimensional lattice for e = 0. 
We show that analytical formula for the exit probabil- 
ity can be derived and the several approaches (among 
them the simple mean field approach) appear to lead 
to the same result. Moreover, received analytical for- 
mula agrees perfectly with Monte Carlo simulations. On 



the one hand, our work is part of a trend of the current 
debate on the form of the exit probability in the one- 
dimensional Sznajd model ^-^5. On the other hand, 
it concerns the much broader problem of nonlinear q- 
voter model. The puzzle, that we deal with here, may 
be contained in a simple question: "Why the mean field 
approach gives the exact formula for the exit probability 
in the one-dimensional nonlinear g-voter model?" 

We have to admit here that our question has been 
strongly inspired by the initial twofold question of Galam 
and Martins "Why the mean field approach gives the 
exact formula for the exit probability in your one- 
dimensional modified Sznajd model or why the Monte 
Carlo simulations give incorrectly a mean field result?" . 
Their recent paper [13] is concluded with "Therefore the 
question is open for future work to settle either there is 
an explanation on why the system studied here exhibits 
a mean field behavior or why different simulations of the 
same system for different sizes all produced incorrectly a 
similar mean field result." We extend their question to 
generalized g- voter model and show by computer simula- 
tions which explanation of a puzzle given in [13| is more 
probable. 



II. NONLINEAR <?- VOTER MODEL IN ONE 
DIMENSION 



We consider a system of L spins 5^ = ±1 located on 
a one dimensional ring. At each elementary time step 
t a panel of q neighboring spins Si, Si-^-i, . . . , Si+q-i is 
picked at random. If all q neighbors are in the same 
state, they influence surrounding spins, if not all spins in 
the q-panel are equal then nothing changes. Two versions 
of the model are considered: 

1. Both sides: The g-panel influences R neighbors 
on the left side and the right side of the panel si- 
multaneously - all spins Si-ii{t + At), Si-ii+i{t + 
At),... Si-i{t + At) and S^+g+i{t + At),Si+q+2it + 
At), . . . Si+q-^-R-l{t + At) take the state of the 



panel, i.e. -^ Si{t). It is easy to notice that for 
R = 1 and q — 2 we deal with the original Sznajd 
model. 

2. Random: The (/-panel influences R neigh- 
bors only on the one, randomly chosen, side 
(left or right) - with probability 1/2 spins 
S.^Rit + At), 5,_H+i(i + Ai), . . . S,-i{t + At) -^ 
Si{t) or ,with the same probability, Si+q+i{t + 
At), Si+g+2it + At),... S^+q+R-i{t + At) -^ S,{t). 
In this case for i? = 1 and q — 2 we deal with the 
modifled version of the Sznajd model, introduced 
by Slanina [11|. For R = 1 and q = 1 we obtain 
original linear voter model 14| . 

After one elementary step time increases by Ai — 1/L. 
Therefore time unit corresponds to one Monte Carlo Step 
(MCS). 

We have introduced both versions to be consistent with 
several other papers that deal with the problem of exit 
probability in the Sznajd model. In some of them "both 
sides" version is used [12[ , whereas others deal with "ran- 
dom" version. However, as we will show, there is no dif- 
ference between the exit probability for both versions and 
therefore each of them can be used. 



III. EXIT PROBABILITY 



scheme [10, llll . Recently it has been shown that even 
much less sophisticated method, simple mean fleld ap- 
proach, leads to the same result [ij]. The question arises 
why the approximate methods give the exact result in 
this case? Some suggestions related to the importance of 
choosing the first pair appeared very recently [l2, Il3| ■ 

In [l2| the Sznajd model of range R (SM(R)) has been 
studied. As usually, at each time step, a pair of near- 
est neighbors Si, S'i+i is chosen at random. If Si = Si+i 
then R neighbors to the left and R neighbors to the right 
change value to Si. This corresponds to the 'both sides' 
version of the nonlinear g- voter model (with q — 2), in- 
troduced in the previous section. Remarkably, in the 
studied case the exit probability E{p) turned out to be 
completely independent of the range of the interaction 
R. Based on this result E{p) was derived in the follow- 
ing way - the exit probability is given by the probability 
that a pair of sites in state -1-1 is chosen before any pair 
of sites in state —1 1121: 



Eip) = p'j2[^-p"-(^-p)"y 



p 



p2 + (1 - p)2 



(2) 



(3) 



Let us consider a finite system with an initial fraction 
p{Q) of randomly distributed spins in the -|-1 state. Two 
scenarios for q- voter model on a one dimensional ring are 
possible: 

• If there is no cluster of size > q, the system is dead- 
locked and no evolution is possible, since only unan- 
imous panel of size q is able to change the state of 
the system. Obviously the number of deadlocks 
grows with q. For q = \ no configuration is dead- 
locked and for q = 2 the only deadlocked configu- 
ration is antiferromagnetic state. 

• If there is at least one cluster of size > q then the 
system will evolve and eventually reach the ferro- 
magnetic state - with probability E{p) the state 
'all spins -|-1' is obtained, whereas with probability 
1 — E{p) the state 'all spins -1' is reached. E{p) is 
called the exit probability and it is one of the most 
important first-passage properties J14l |. 

In 2008 the exit probability for the one-dimensional Sz- 
najd model (which corresponds to nonlinear voter model 
with q = 2) has been calculated analytically [ifl, llil : 



E{p)^ 



P 



p2 + (1 - p)2 



(1) 



This result agrees perfectly with Monte Carlo simula- 
tions, which is quite puzzling since calculations were not 
exact but based on Kirkwood-approximation decoupling 



Another idea was presented in . 13| - the quasi- 
deterministic procedure in which the only random step 
is the selection of the initial pair, and then the process 
is deterministic. These types of considerations have led 
again to the same formula for the exit probability ([1]). 

In this paper we show that mean field approach gives 
the formula that agrees perfectly with Monte Carlo sim- 
ulations for arbitrary value of q and we also examine pos- 
sible explanations why approximate methods give the ex- 
act result for a one-dimensional nonlinear g- voter model. 

We start with showing by computer simulations that 
exit probability for both version of the model ( "random" 
and "both sides") are identical. In Fig. [T] we present 
results for the lattice size L = 100 and two values of 
q. Analogous results have been obtained earlier for the 
Sznajd model (i.e. q = 2) [il|. Given that both versions 
of the model give identical results, we will concentrate in 
the later work on "random" version. However, we have 
checked all results presented in this paper for both cases. 

Now we are ready to derive analytical formula for q- 
voter model. Let us denote by < N^{t) > the average 
value of the number of g-panels with all spins -|-1, and 
by by < N~{t) > the average value of the number of 
g-panels with all spins +1. Obviously: 



<N+it)>^pit)^L 
<N^{t)>={l~p{t))'^L, 



(4) 




both sides, q=3 
random, q=3 
both sides, q=5 
random, q=5 
-analytical formula 




FIG. 1: Exit probability in the case of the "random" and 
"both sides" nonhnear q-voter model on a one-dimensional 
ring of length L = 100. Averaging was done over 10* samples. 
Results are identical for both versions of the model. 



FIG. 2: Exit probability for nonlinear q-voter model on a one- 
dimensional ring of length L = 100. Analytical formula agrees 
with the Monte Carlo results for any value of q. Averaging 
was done over 10* samples. 



where p{t) is the fraction of spins +1 at time t. Let us 
now introduce the quantity; 



mq{t) 
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^; 


{t)>-< N-q 


-{t)> 


< 


N,^ 


{t)> + < Ng 


-{t)> 
- 1. 



p1{t) + {l-p)1it) 



(5) 



which for q = 1 (linear voter model) is simply average 
magnetization: 



mi(t) 



<N+{t) > - <N^it) > 

L 
2p{t) - 1. 



(6) 



It is known that for linear voter model average magne- 
tization is constant, i.e. mi{t) — r7T,i(0) [ij]. Knowing 
this, it is very easy to derive exit probability: 

m,(oo) = Eip)-{l-E{p))=2Eip)-l, (7) 

where we use notation p{0) = p. Because 

TOi((X)) = mi(0) = 2p- 1, (8) 

we obtain that: 

E{p) = p. (9) 

Now we could generalize this reasoning assuming that: 

mg{oo) = mg{0). (10) 

With such an assumption we obtain result for arbitrary 
value of q: 



E{p) 



P" 
p<i + {l-p)i- 



(11) 



Indeed substituting to (|lll) value of g = 1 we obtain 
known result E{p) — p and for q = 2 we obtain formula 
((TJ derived independently in 4 papers |10l - [l3| . Of course, 
since we did not show any proof that our assumption 
PH)) is valid, the formula pT|) can be treated as a guess. 
Let us start with checking validity of the formula pT|) by 
performing Monte Carlo simulations. In Fig[2]we present 
both, analytical and Monte Carlo, results for several val- 
ues of q. As seen there is a perfect agreement though the 
number of averaging is not very large (10* samples). 

Now we will follow the reasoning presented in [I4I . We 
check how the exit probability depends on the range of 
interactions R. Up till now we focused only on i? = 1 , but 
we have provided simulations also for R = 2,3,4,5 and 
10. Results for i? = 1 and R = 5 are presented in Fig. 
131 Analogous results were obtained also for other values 
of R. It occurred that exit probability for the g-voter 
model is identical for any value of R. The same results 
has been obtained for the Sznajd model (i.e. q = 2) in 
[1^. Remarkably, also the size of the system does not 
influence the exit probability (see Fig. [3]). The same 
results have been obtained earlier for the Sznajd model 
[111 Il2| . This is quite intriguing, since in the case of a 
g- voter model on a complete graph the size of the system 
influences the exit probability -E(/o)[M|- Using the fact 
that the range of interactions does not change the exit 
probability we can derive analytical formula for E{p) in 
the same way as in [i2|. For R> {L — q)/2 (in the case 
of "both sides" ) ot R> L — q (for "random" ) the system 
is fully ordered after first unanimous g-panel is chosen. 
Therefore the exit probability is equal the probability of 
choosing g-panel of 'up' spins: 



E{p) 



P" 

p9 + (l-p)9' 



(12) 
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FIG. 3; Exit probability for nonlinear (q=2)-voter model on 
a one-dimensional ring of length L. Exit probability changes 
neither with the lattice size L nor with the range of interac- 
tions. Analogous results has been obtained in [l(]| - [l3 |. Anal- 
ogous results are valid for any value of q. Averaging was done 
over 10'' samples. 



FIG. 4: Difference between two exit probabilities E^{p) and 
E~ (p) for the lattice size L — 100 and the range of interaction 
R = 1. In the case of E^{p) a cluster of q 'up' spins has been 
chosen in the first step of simulation, whereas in the case 
of E~ (p) a cluster of q 'down' spins has been chosen in the 
first step. It is seen that the first choice influence the exit 
probability. Averaging was done over 10^ samples. 



This type of reasoning coincides, in a sense, the idea 
presented in the work of Galam and Martins [Tj] . They 
have proposed simple quasi-deterministic procedure that 
drives system to the absorbing ferromagnetic state with 
the same exit probability as in the Sznajd model. It their 
procedure the only probabilistic step is the choice of the 
first pair and further evolution is entirely deterministic. 
To see whether the first choice actually affects the proba- 
bility of the final state we performed two types of simula- 
tions. In the first step we were picking at random g-panel 
of 'up' spins and then the system was evolving under 
standard procedure or in the first step we were picking 
at random g-panel of 'down' spins. We have denoted the 
exit probability in the first case by E^{p), whereas in the 
second case E~^{p). If the first choice does not infiuence 
the exit probability we should obtain E^{p) = E^{p). 
We have plotted in Fig. |3]the difference E^{p) — E^{p) 
for several values of q. As seen the difference grows with 
q, which is quite understandable. For large value of q 
this is difficult to find unanimous g-panel, especially for 
/9 — 7- 0.5 (where only small clusters are present in the 
system). Probability of finding g-panel of 'up' spins is 
equal p'^, whereas probability of g-panel of 'down' spins 
is equal (1 — p)'^ . This explains shapes of the curves in 
Fig. 21 However, for g = 2 importance of first choice is 
not too high and therefore quasi-deterministic procedure 
proposed by Galam and Martins [Ij] cannot explain the 
analytic formula for E{p). Also the results concerning 
the importance of the interaction's range R should be 
treated rather statistically. The fact that R does not in- 
fiuence the exit probability does not mean that it is not 
important at all. In other case the first choice would de- 
termine the final state of the system, which is obviously 
not true (see Fig. |3]and[5]). 




FIG. 5: Difference between two exit probabilities E (p) and 
E^ (p) for several lattice sizes L and the range of interaction 
R = 1. In the case of E^{p) a cluster of g = 2 'up' spins 
has been chosen in the first step of simulation, whereas in the 
case of E~ (p) a cluster of g = 2 'down' spins has been chosen 
in the first step. Results scale trivially with the system size 
L. Averaging was done over 10^ samples. 



Although the choice of the first pair does not determine 
the result in 100% but the first steps of simulation are 
actually the most important. To see this let us present 
probability p{t) of choosing g-panel of 'up' spins as a 
function of initial probability p for several values of t 
(see Fig. [6]). As we see, already after several Monte Carlo 
steps p{t) concides with E{p). Why the latter steps does 
not change the form of E{p)7 

To understand this, we should mention here that the 
formula (llip is valid only for random initial conditions. 




FIG. 6: Fraction of 'up' spins after time t measured in the 
Monte Carlo Steps. 



If we start from two clusters - the first cluster of spins 
+1 and length pL and the second one of spins —1 and 
length (I — p)L we obtain the final state of all spins +1 
with probability p. This is easy to understand because 
the probability in such a case of choosing g-panel of 'up' 
spins is equal p (at least for the infinite system size L). 
The same result has been obtained for the Sznajd model 
in more sophisticated way using Kirkwood approach [11| . 



Therefore, as soon as the system will order to several 
domains then E(p) will no longer continue to change. 



IV. SUMMARY 

A recent work [l3| did point out to the evoked puz- 
zle emphasizing the mean field character of the formula 
for the exit probability in the Sznajd model obtained in 
[lOl lll| . It has been shown that the continuous shape of 
the exit probability is a direct outcome of a mean field 
treatment [l^] ■ Two possible explanations has been given 
- most likely finite size effects in the simulations or as 
a alternative irrelevance of fluctuations in the system. 
In this paper we have extend the puzzle to the nonlin- 
ear g-voter model on a one-dimensional lattice. We have 
proposed analytical formula for the exit probability and 
checked it validity by Monte Carlo simulations for sev- 
eral values of q. In particular, we have shown that neither 
the range of interactions nor the size of the system influ- 
ence the exit probability. Moreover, we have checked the 
importance of the initial evolution of the system, as sug- 
gested for the Sznajd model in [Tj, [lj|, to understand 
proposed analytical formula. It should be stressed that 
results presented here contain a well-known cases of lin- 
ear voter and Sznajd model. 
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